Rotational damping in a multi-j shell particles-rotor model 
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Abstract 

The damping of collective rotational motion is investigated by means of particles-rotor model in 
which the angular momentum coupling is treated exactly and the valence nucleons are in a multi-j 
shell mean-held. It is found that the onset energy of rotational damping is around 1.1 MeV above 
yrast line, and the number of states which form rotational band structure is thus limited. The 
number of calculated rotational bands around 30 at a given angular momentum agrees qualitatively 
with experimental data. The onset of rotational damping takes place gradually as a function of 
excitation energy. It is shown that the pairing correlation between valence nucleons has a significant 
effect on the appearance of rotational damping. 
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I. INTRODUCTION 



The experimentally observed rotational bands often lie in the region near yrast line and 
are described as particle-hole excitations in the mean-field. At higher excitation energy above 
the yrast line, it does not necessarily form rotational band structure due to the damping of 
collective rotational motion . When the rotational damping takes place, E2 transition 
from an excited state spreads out over many final states. The gamma-rays which are emitted 
from the above excited region can not be distinguished as discrete peaks, thus forming a 
quasi- continuum spectra. Experimentally rotational damping has been studied through 



the analysis of quasi-continuum spectra 







. Recent experimental progress in high 



precision three-dimensional gamma-ray correlation measurements makes it possible to study 
various features of collective rotational motion in the regions of discrete rotational bands and 



damped rotational motion 8 
fluctuation analysis method 
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151 ] . In particular, the newly developed 



8J, |9j] has presented the number of rotational bands existing in 
a rare-earth nucleus is only about 30 at a given angular momentum, thus confirming the 
occurrence of rotation damping. 

Early theoretical studies on rotational damping showed that with the in- 

crease of level density the off-diagonal residual interaction becomes effective to cause mixing 
of many-particle many-hole configurations in the rotating mean-field. Since different con- 
figurations respond differently to the Coriolis force, the configuration mixing results in a 
dispersion of rotational frequency within each energy eigenst ate, implying a damping of the 
collective rotational motion. However, in these works (3, 0, Q, 0|, the assumption that 
configuration mixing is described by the general statistical theory of random matrices has 
been used to treat the E2 strength function associated with the damped rotational motion. 

Recent studies on the microscopic structure and mechanism of rotational damping have 
been done extensively for normally deformed and superdeformed nuclei in the cranked Nils- 
son mean-field combined two-body residual interaction 

these discussions, the excited rotational bands are described as intrinsic many-particle many- 
hole excitations in the cranking Nilsson mean-field. In order to obtain rotational damping, 
the bands are mixed by two-body residual interaction. Furthermore, the cranking model is 
semi-classical and the angular momentum is not a good quantum number. 

One can also use the particles-rotor model (PRM) to obtain the band structure as well 
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as rotational damping and to study the coupling between single-particle degree of freedom 
and collective rotational motion. Moreover as noted in Ref. [l9|, PRM is more appropriate 
on angular momentum coupling in wave functions and E2 transition matrix elements, since 
angular momentum is not a conserved quantity in cranking model. With many particles 
in a single- j shell, a systematic comparison of cranking model and PRM has been done 
and the pair-correlation transition in rotating nuclei has been investigated in Ref. 26^ . 
In Ref. [27], PRM has been used to examine the quality of the tilted axis cranking and 
its interpretation. In particular, by correctly treating the angular momentum coupling in 
triaxial PRM, a new phenomenon-chiral doublet bands has been predicted [3] and verified 
later by experiment. Therefore it will be interesting to study the rotational damping by 
PRM. The nuclear chaotic behavior and its connection with rotational damping have been 
addressed in cranking model and PRM with single- j shell 



30, 



3 it ] . The purpose of present 



work is to study the rotational damping by PRM, in which the angular momentum is a good 
quantum number and the pairing correlation will be included explicitly. The mechanism of 
rotational damping and E2 transition property from the region of discrete rotational band 
into the regime of damped rotational motion will be investigated in a multi-j shell PRM. 



II. FORMULATION 



A nucleus can be visualized as a rotor coupled with a few valence nucleons, which move 
more or less independently in the deformed potential of core composed by the rest of the 
nucleons, i.e., the particles-rotor model (PRM). The Hamiltonian of PRM is expressed as 
the sum of an intrinsic part and a collective part 

Hprm = H intv + H co \ h (1) 

where Hi ntr describes microscopically the motion of valence particles near Fermi level and 
H co \\ is the collective rotation of core. The intrinsic Hamiltonian is taken as 

#intr = 0'l m l| ~ ^ K -V^ Y 20\j2m 2 )a] imi a j2m2 - G ^ a ]'m' a ]'fh' a j^ a jm, (2) 

jiTnij2TTi2 j'm'jm 

where aj m and the one-particle creation and annihilation operators and a residual 

pairing interaction has been included. In order to describe many-particle many-hole excita- 
tions associated with rotational damping, we include some different- j shells in the deformed 
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mean-field of valence nucleons. The single-particle energy in the deformed mean-field is 



written as 129 
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j\m\j%m2 j m 

where Rj is a parameter to indicate the relative energy between different- j shells. The 



deformation parameter k is related to the quadrupole deformation (3 through 
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« ~ 0.16/iu;o(iV + 3/2)/3, (4) 

where ftwo is harmonic oscillator frequency of the deformed potential and iV the quantum 
number of the major shell. For example, in the case of ii3/ 2 shell, k = 2.5 MeV approximately 
corresponds to (3 ~ 0.3, and h§/ 2 shell corresponds to k — 2.2 MeV. The two-body correlation 
between valence particles is taken as pairing interaction with strength parameter G. Since 
we solve two-body interaction exactly, it will contribute to both particle-particle channel 
and particle-hole channel. 

The spin / of nucleus is the sum of angular momentum R of core and J, the sum of 
angular momentum of valence nucleons. The collective Hamiltonian can be expressed as 

R 2 + R 2 _ I 2 -I 2 J 2 -J 2 J + J_ + /-J + 
co11 2J 2J + 2J 2J 1 } 

where H TOt = (I 2 — J|)/2J7" describes the nuclear collective rotation. The recoil term H rec = 
(J 2 — J 2 )/2J~ acts only on valence nucleons and contains one-body and two-body terms 
if there is more than one particle outside the core. Coriolis interaction H cor = — (J + J_ + 
I-J + )/2J couples the collective rotation and the single particle motion. 
The eigenf unctions of particles-rotor Hamiltonian can be expanded as 

*?M = J2 C KrfMK (6) 
K 

with the mixing coefficient C^, where the symmetrized basis is given by 



T> I MK (Vt) is the usual rotation matrix and (j>^a ^ e ^-body anti-symmetrized wave func- 
tion of valence particles. 
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The stretched E2 transition probability from an initial nuclear state a at spin / + 2 to a 
final state a' at I is 

B(E2, aI + 2^ a' I) = ^Q 2 M 2 aI+2tCefI , (8) 

where the amplitude 

M aI+2 , a ,j = J2 c k{I + 2)C&(I)(IK'20\I + 2K)8 KK , (9) 

KK' 

We neglect the minor contribution of non-collective E 2 transition from valence particles 
since the electric quadrupole moment Qq is much larger than that from valence particles. The 
CG coefficient (IK20\I + 2K) represents the coupling of angular momenta in the intrinsic 
frame. The normalized E2 transition probability S a i + 2, a 'i is defined as 

S a I+2,a'I = M ll+2,a'n 

which satisfies J2 a ' Sai+2,a'i = 1 from an initial state a to many final states a'. 



III. RESULTS AND DISCUSSIONS 



dls. These pc 



In present PRM, we consider eight valence particles in a multi-j shell coupled with a de- 
formed core. The calculation has been performed with the parameters: the pairing strength 
G = 0.45 MeV, the moment of inertia ,7=76.0 h 2 MeV~\ k = 2.5 MeV for i 13/2 intruder 
orbit and k = 2.2 MeV for other shells. These parameters have been used in earlier studies 
and are considered to be reasonable 

In order to imitate the level scheme of Nilsson diagram |35|, we include the shells h 9 / 2 , 
H3/2, P3/2, f$/2 and pi/2 in the deformed mean-field. The relative energies Rj between 
different- j shells are determined according to spherical Nilsson level diagram. Table H] gives 
the relative energy Rj between shell h 9 / 2 and other shells H3/2, P3/2, /s/2 and P1/2, e.g. it 
has been taken as 0.1 Huq between shells hg/2 and ii3/2- The choice of these parameters for 
the present single-particle levels is similar as the neutron levels in nucleus 168 Yb. Of course, 
the valence proton excitations should be included as well to fully understand nucleus 168 Yb. 
However, if both proton and neutron excitations are taken into account, the Hamiltonian 
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281 ]. only one proton and one 



matrix will be too big to be diagonalized exactly. In Refs. 
neutron are coupled with the rotor, and the many-particle and many-hole excitations are 
thus missing. 



The level scheme thus obtained is shown in Fig. ^ in which the spherical shells and their 
splitting due to deformation are presented. On the right the third component rrij of angular 
momentum is given. 

After defining the single-particle basis, the many-particle many-hole excitations within 
a given excitation energy form the shell model basis states. The configuration space is 
spanned by various excitations of eight valence particles in the deformed mean-field. The 
off-diagonal Hamiltonian from pairing correlation, recoil term and coriolis interaction will 
mix states not only in the same-j shell, but also in different-j shells 0, EiZ|- They are 
responsible for the many-particle many-hole excitations. In the model diagonalization for 
each spin / = — 60 h, the configuration truncation is done with the excitation energy ~ 0.6 
fkjQ, which corresponds to the lowest 2000 many-particle many- hole configurations. Most of 
these configurations has lp — lh to 4p — Ah characters and few of them is with 5p — 5h to 
6p — 6h characters. 

Fig. 121 displays the calculated nuclear states with small horizontal bars. Strong E2 tran- 
sitions satisfying condition S a i + 2 ia 'i > 1/V2 = 0.707 are plotted with solid lines. Weaker 
transitions defined by condition 0.5 < S a i + 2,a'i < 0.707 are presented with dashed lines. 
One may observe that most of the strong E2 transitions lies in the region near yrast line, 
where the rotational band structure is identified. At higher excitation energy the transitions 
become much weaker and E2 transition from an initial state may spread out over many final 
states, implying the disappearance of rotational band structure. 

In order to understand the nature of nuclear state, a quantity, occupation number of 
single-particle basis [Ijrrij] in nuclear state a, is defined as 

a(ljnij) 

satisfying condition that the sum of various occupancy on single-particle basis is equal to 
the number of valence particles 

J2 P ^=N. (12) 

Ijnij 

Here we use quantum numbers [Ijmj] to denote the single particle basis states in Fig. 
Fig. El gives the occupation number of particle for various single-particle basis [Ijrrij] as a 
function of angular momentum for (a) ground state band; (b) a strong transition band at 
excitation energy U ~ 0.8 MeV; (c) a weak transition band at excitation energy U ~ 1.6 
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MeV. The lines with filled upper triangle, cross, filled circle, open square, open circle, filled 
squire and open upper triangle represent the single-particle basis states [5||], [6y|], [5|§], 
[6y|]> [^yIIj [5ff] an d [6y|] separately. Since the occupancy on other higher single 
particle basis is not important for the low excited states, we use the lines with star to 
represent the occupancy on those single particle basis. One may observe for the bands 
satisfying strong transition condition as shown in Fig. El (a) and (b), the occupation number 
of particle is nearly independent of angular momentum. With the increase of excitation 
energy, the dependence of occupancy on spin becomes more sensitive, as shown for the band 
with weak transition in Fig. El (c). Since the occupancy characterizes the wave function 
property of nuclear state, these figures clearly indicate how the rotational band structure 
starts to disappear gradually due to the change of wave function in nuclear states. 

In the experiments, strong E2 transitions are observed as discrete peaks in the gamma- 
ray spectra and the rest of transitions shows up as quasi-continuum spectra which contain 
transitions summed over many final states. For such a situation, it is useful to represent the 
E2 transition property by means of strength distribution function. The strength function 
for a state a at spin / + 2 is given by [19] 

S a ,I+2(E~ / ) = S a I+2,a'l5(E~ / — E a j +2 + E a 'j). (13) 

a' 

The fragmentation of E2 strength function is the rotational damping phenomenon. In order 
to quantify the onset of rotational damping, the branch number is defined as 

^branch {al + 2) = S^+^j) _1 . (14) 

a' 

It counts effectively the branch number of E2 transitions from an initial state a at / + 2 
to the final states a', which are allowed by the selection rules of the gamma-ray radiation. 
For a case where a given state decays to only one final state at the lower spin, branch is 
equal to one. If there are two possible transitions with equal probability, the branch number 
is equal to two. In other words, nt, ra nch < 2 implies that the transition from a given state 
mainly decays to one final state, where the discrete rotational band structure is identified. 
In contrast, branch > 2 states that the transitions spread over two or more final states, 
where the rotational damping takes place. 

In order to indicate the E2 transition property more precisely, the strength function 
S a j(E^) is shown in Fig. 0] for the lowest 9 states with I n = 30 + . The branch number 
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^branch and excitation energy U are put for each state. The E2 transition associated with 
the first 30 + state exhausts most of the total strength at gamma-ray energy 0.92 MeV. The 
30^, 30^, 30J and 30^~ states show essentially the same E2 distributions except for slight 
difference in the gamma-ray energy. The E2 decays from the 30g~, 30^", 30^ and 30g~ states 
display a different E2 strength distribution, being fragmented over several transitions, each 
of which carries a rather weak strength. Fig. EJdisplays the quantity S , Cfj j-(£' 7 ) associated with 
the states 30g 5 and 30g 6 lying at excitation energy U ~ 2.3 MeV, and the states 30g 7 and 30<jg 
lying at U ~ 3.1 MeV. The E2 strength distribution at U ~ 2.3 MeV has about 16 branches, 
while the number of branches becomes around 42 at U ~ 3.1 MeV. It indicates that the 
degree of the mixing between many-particle many-hole configurations becomes stronger as 
excitation energy increases, which is the basic mechanism to form the damping of collective 
rotational motion. 

As the excitation energy increases, the rotational band structure gradually disappears 
and rotational damping takes place. Branching number is a key quantity to measure where 
rotational damping takes place and the degree of configuration mixing. The dependence 
of branching number on excitation energy is shown in Fig. IHlfor spins (a) I n = 20 + , (b) 
I w = 30 + and (c) I 77 = 40 + . The histogram gives the average branch number within the 
energy bins. One may observe that branching number increases with excitation energy. 
Using the criterion branch = 2 for the onset of rotational damping, the onset energy is 
predicted to be around excitation energy U ~ 1.1 MeV above yrast line. 

It should be noted, although the onset energy defined by condition branch = 2 tells 
approximately where the rotational damping takes place, the transition from the region of 
rotational bands to rotational damping is not very sharply at the onset energy, but develops 
gradually as the excitation energy increases. As shown in Fig. |21 there exist some rotational 
bands in the region of rotational damping, where the excitation energy is higher than the 
onset energy. These band structures are surrounded by states which do not have strong 
transitions. It indicates that the rotational band structure partly remains even in the region 
of rotational damping. Such feature of rotational bands is also displayed in Fig. |HJ where 
there exist states whose branching number is smaller than 2 at higher excitation energy. 

Due to the onset of rotational damping, the number of rotational bands existing in a 
given nucleus is thus limited and gives a quantitative measure of rotational damping. Here, 
the number of rotational bands corresponds to the experimental effective number of paths 
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which can be obtained from x E 7 spectrum by the fluctuation analysis method 
Theoretically the number of rotational bands is defined by the number of states with strong 
E2 transition, which satisfies condition S a i + 2,a'i > 0.707 or n brancn < 2 with at least two 
consecutive steps 1 + 2^1^1 — 2ofE2 decays [3]. Since the configuration of valence 
nucleons occupancy in present model corresponds to rare-earth nucleus 168 Yb considered as 
a core 160 Yb plus eight valence neutrons and in our calculation the moment of inertia is 
taken as the value of rare-earth nucleus, we will make a qualitative comparison between our 
calculated number of rotational bands and the experimental data for nucleus 168 Yb, though 
the model space in present calculation is not really realistic for nucleus 168 Yb. For this 
nucleus there exists experimental data from the analysis of quasi-continuum gamma-spectra 
as well as data from discrete spectra identifying the rotational bands up to spin / ~ 40 
j^ls^. Fig. [7| shows the calculated number of rotational bands as well as the experimental 
effective number of paths P| for nucleus 168 Yb. The solid line represents the result with the 
criterion S a i + 2, a 'i > 0.707 while the dashed line is calculated with condition nbranch < 2. 
The horizontal axis denotes the average gamma-ray energy = (E 7l + E j2 )/2, where E ll 
is the transition energy for / + 2 — > I and E 12 for transition I — > / — 2. It is clear that 
the two conditions give essentially the same number of rotational bands around 30, and the 
theoretical calculation agrees well with the experimental result in all the gamma-ray energy 
range. 

In order to study the effect of pairing correlation between valence nucleons on rotational 
damping, Fig. [Hldisplays the calculated number of rotational bands with and without pairing 
correlation, together with experimental result. The solid line represents the result without 
pairing interaction between valence nucleons, while the dashed line is with the standard 
pairing G = 0.45 MeV. The criterion S a i + 2,a'i > 0.707 has been used to obtain the number 
of rotational bands. It is seen that the calculated number of bands becomes much larger when 
the pairing is turned off. In other words, the rotational band structure has strengthened and 
rotational damping has weakened when pairing G = 0.0 MeV. This indicates that pairing 
correlation has a significant effect on the damping of rotational motion. 

To understand the nature of nuclear state without pairing correlation, Fig. |§] gives the 
occupation number of particle of the same bands as in Fig. |3J but with pairing strength G = 
0.0 MeV. One may observe for these bands the dependence of occupancy on spin becomes 
weaker in the case of without pairing correlation. Especially in low angular momentum, 
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there is much difference between the occupancy with and without pairing. These results 
are consistent with the calculation in Fig. |H1 where the number of rotational bands becomes 
larger when paring G = 0.0 MeV. It indicates that pairing correlation has an important 
effect on the nature of nuclear states especially in low angular momentum. 

It should be mentioned that pairing correlation between valence nucleons has contribu- 
tions to both diagonal and off-diagonal Hamiltonians. Diagonal Hamiltonian characterizing 
the property of mean-field favors the rotational band structure and retards the rotational 
damping, whereas the off-diagonal Hamiltonian characterizing the quantum fluctuations 
coming from residual interaction causes the mixing of many-particle many-hole configu- 
rations and prefers the rotational damping. Therefore, the delicate balance between the 
competition of diagonal and off-diagonal components of pairing correlation has determined 
the appearance of rotational damping. It is concluded that the off-diagonal components of 
pairing correlation (two-body residual interaction) play an important role in the appearance 
of rotational damping. 



IV. CONCLUSIONS 



The damping of collective rotational motion is discussed in a multi-j shell particles-rotor 
model, in which the angular momentum is strictly conserved and the pairing correlation is 
included explicitly. It is found that the rotational damping takes place at about 1.1 MeV 
above yrast line, and the number of states which form rotational band structure is thus lim- 
ited. The onset energy in present calculation is similar as the theoretical prediction U ~ 0.8 
MeV in cranked Nilsson mean-field combined two-body residual interaction. The calculated 
number of rotational bands around 30 is in a qualitative agreement with experimental result 
in the gamma-ray energy range. The onset of rotational damping takes place quite gradu- 
ally as a function of excitation energy. Even in the region of rotational damping, there still 
remains part of discrete rotational band structure. Our calculation clearly indicates that the 
pairing correlation has an important effect on the nature of nuclear states, especially in low 
angular momentum. It is found that the calculated number of bands becomes much larger in 
the case of pairing strength G = 0.0 MeV. The pairing correlation between valence particles 
has significant effect on the appearance of rotational damping. It is noted that the pairing 
correlation and the exact treatment on angular momentum coupling in PRM play an impor- 
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tant role to provide a description on rotational damping. However in the cranked Nilsson 
mean-field combined two-body residual interaction, the rotational damping is attributed to 
the high-nnltipole component of two-body residna, tntetaction Q. Considering that the 
angular momentum in Ref. [19| is not treated properly and the valence nucleons states are 
schematic in present model, more works along these lines are necessary to understand better 
the mechanism of rotational damping. 
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TABLE I: The relative energy Rj between shell hg/ 2 and other shells £13/2, P3/2) fh/2 an d P1/2 i n 
unit of fuvQ. The level scheme is shown on the left in Fig. ^ 
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FIG. 1: Single-particle level scheme in the mean-field of valence particles. The spherical shells on 
the left split up to the structure on the right due to the deformation with the third component rrij 
of angular momentum. 



14 





10 20 30 40 50 60 

Spin (h) 

FIG. 2: The calculated nuclear states are plotted with small horizontal bars. A reference energy 
I{I+1)/2J with ^7=76.0 7i 2 MeV _1 is subtracted. Solid lines connecting the energy levels represent 
the strong E2 transitions which have the normalized strength S a i+2,a>i larger than 0.707. Dashed 
lines are the weaker E2 transitions with normalized strength between 0.5 and 0.707. 
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(a) ground state band 




















20 30 40 50 60 

Spin (ft) 

(b) a strong transition band at J/-0.8 MeV 
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(c) a weak transition band at U-l.6 MeV 
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FIG. 3: Occupation number of particle for various single-particle basis [Ijrrij] as a function of spin 
for (a) ground state band; (b) a strong transition band at excitation energy U ~ 0.8 MeV; (c) 
a weak transition band at excitation energy U ~ 1.6 MeV. The lines with filled upper triangle, 
cross, filled circle, open square, open circle, filled squire and open upper triangle represent the 
single-particle basis states [5§|], [6^±], [5§§], [6^|], [6^|], [5§§] and [6^|] separately. The 
lines with star denote the occupancy on other single particle basis. 
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FIG. 4: The E2 strength distribution S^E^) from the lowest 9 states with I n = 30+ . The branch 
number ribranch and excitation energy U are put for each state. 



0.20 



0.15 



0.10 



0.05 



w 0.08 



0.06 
0.04 
0.02 



0.00 



35-th 

U=2.28 MeV 

nbranch=15.21 



97-th 

U=3.08 MeV 

n b ranch=42.34 



36-th 

U=2.31 MeV 

nbranch=17.08 



98-th 

U=3.11 MeV 

ntjranch=41 .48 



0.20 
0.15 
0.10 
0.05 
0.08 
0.06 
0.04 
0.02 
0.00 



0.6 0.8 1.0 1.2 1.4 1.6 0.8 1.0 1.2 1.4 1.6 

£ 7 (MeV) £ 7 (MeV) 



FIG. 5: The strength distribution S a j(E 1 ) for the stretched E2 decays from the 35th and 36th, 
97th and 98th excited states with r = 30+. 
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FIG. 6: The branching number n branch for (a) F = 20+, (b) J" = 30+ and (c) F = 40+ as a 
function of excitation energy U. The histogram gives the average branch number within the energy 
bins. The horizontal line shows nb ra nch=2 used to define the onset of rotational damping. 
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FIG. 7: The calculated number of bands together with the experimental effective number of paths 

Q 

for nucleus 168 Yb _9] as a function of the average transition gamma-ray energy. The solid line is 
calculated with the criterion S a j + 2,a'i > 0.707 while the dashed line represents the result with 
condition ribranch < 2. 
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FIG. 8: The calculated number of bands with and without pairing correlation together with the 
experimental result as a function of the average gamma-ray energy. The solid line is the result 
with pairing strength G = 0.0 MeV, while the dashed line is with the standard pairing G = 0.45 
MeV. The criterion S a i + 2,a>i > 0.707 has been used to obtain the number of rotational bands. 
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FIG. 9: The same as in Fig. |31 but without pairing correlation. 
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